The hagfish genome and the evolution of vertebrates.

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Gene expression in the liver of the hagfish (Eptatretus burgeri) belonging to the Cyclostomata is ancestral to that of mammals.

Although the liver of the hagfish, an earliest diverged lineage among vertebrates, has a histological architecture similar to that of mammals, its gene expression has not been explored yet. The present study was undertaken to comparatively characterize gene expression in the liver of the hagfish with that of the mouse, using in situ hybridization technique. Expression of alb (albumin) was detectable in all hepatocytes of the hagfish liver, but was negative in intrahepatic bile ducts. Their expression in abundant periportal ductules was weak. The expression pattern basically resembled that in mammalian livers, indicating that the differential expression of hepatocyte markers in hepatocytes and biliary cells may have been acquired in ancestral vertebrates. alb expression was almost homogeneous in the hagfish liver, whereas that in the mouse liver lobule was zonal. The glul (glutamate-ammonia ligase) expression was also homogeneously detectable in hepatocytes without zonation, and weakly so in biliary cells of the hagfish, which contrasted with its restricted pericentral expression in mouse livers. These findings indicated that the hagfish liver did not have mammalian-type zonation. Whereas tetrapods had Hnf (hepatocyte nuclear factor) 1a and Hnf1b genes encoding the transcription factors, the hagfish had a single gene of their orthologue hnf1. Although HNF1α and HNF1ß were immunohistochemically detected in hepatocytes and biliary cells of the mouse, respectively, hnf1 was expressed in both hepatocytes and biliary cells of the hagfish. These data indicate that gene expression of hnf1 in the hagfish liver may be ancestral with that of alb and glul during vertebrate evolution.

Bioinformatic and fine-scale chromosomal mapping reveal the nature and evolution of eliminated chromosomes in the Japanese hagfish, Eptatretus burgeri, through analysis of repetitive DNA families.

In the Japanese hagfish, Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n = 52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n = 36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome of E. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1-6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes of E. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7-10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1-10) by fluorescence in situ hybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2-10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

High-quality total RNA extraction from early-stage lamprey embryos.

High-purity total RNA extraction from animal embryos is essential for transcriptome analyses. lampreys, together with hagfish, are the only extant jawless vertebrates or cyclostomes and are thus key organisms for EvoDevo studies. However, extracting uncontaminated RNA from early-stage embryos remains challenging. RNA does not bind to the silica membrane in filter-based extractions, significantly reducing yields; and ethanol/isopropanol precipitation methods lead to contaminants, bringing down the optical density (OD) 260/280 ratio. The RNA extraction protocol was modified using precentrifugation and adding salts before isopropanol precipitation. This modification significantly increased RNA yield, removed contaminants and improved RNA integrity. Egg membrane sources were suspected to cause RNA purification problems because low-quality extraction does not occur in posthatching embryos.

Phylogenetic and functional properties of hagfish neurohypophysial hormone receptors distinct from their jawed vertebrate counterparts.

Vertebrate neurohypophysial hormones, i.e., vasopressin- and oxytocin-family peptides, exert versatile physiological actions via distinct G protein-coupled receptors. The neurohypophysial hormone receptor (NHR) family was classically categorized into four subtypes (V1aR, V1bR, V2R and OTR), while recent studies have identified seven subtypes (V1aR, V1bR, V2aR, V2bR, V2cR, V2dR and OTR; V2aR corresponds to the conventional V2R). The vertebrate NHR family were diversified via multiple gene duplication events at different scales. Despite intensive research effort in non-osteichthyes vertebrates such as cartilaginous fish and lamprey, the molecular phylogeny of the NHR family has not been fully understood. In the present study, we focused on the inshore hagfish (Eptatretus burgeri), another group of cyclostomes, and Arctic lamprey (Lethenteron camtschaticum) for comparison. Two putative NHR homologs, which were previously identified only in silico, were cloned from the hagfish and designated as ebV1R and ebV2R. In vitro, ebV1R, as well as two out of five Arctic lamprey NHRs, increased intracellular Ca2+ in response to exogenous neurohypophysial hormones. None of the examined cyclostome NHRs altered intracellular cAMP levels. Transcripts of ebV1R were detected in multiple tissues including the brain and gill, with intense hybridization signals in the hypothalamus and adenohypophysis, while ebV2R was predominantly expressed in the systemic heart. Similarly, Arctic lamprey NHRs showed distinct expression patterns, underscoring the multifunctionality of VT in the cyclostomes as in the gnathostomes. These results and exhaustive gene synteny comparisons provide new insights into the molecular and functional evolution of the neurohypophysial hormone system in vertebrates.

A new species of six-gilled hagfish (Myxinidae: Eptatretus) from the Lakshadweep Sea.

A new species of hagfish, Eptatretus wadgensis sp. nov., is described from the Wadge Bank, Lakshadweep Sea, India, obtained from a depth of ~250300 m through deep-sea trawling. It is diagnosed by having six pairs of gill pouches and gill apertures, 3/3 multicusp teeth, total slime pores 6769, six branchial slime pores, and ventral aorta bifurcating at the 4th or between 4th and 5th gill pouch. The new species has significant morphological differences in total dental cusps, total slime pores, body proportions and the absence of the nasal-sinus papilla when compared to congeners and formed a distinct clade in phylogenetic reconstruction and a genetic distance of 3.414.00% when comparing K2P parameters with the nearest species. A key to the Eptatretus species of the Indian Ocean is provided.

A glimpse into the past: phylogenesis and protein domain analysis of the group XIV of C-type lectins in vertebrates.

BACKGROUND: The group XIV of C-type lectin domain-containing proteins (CTLDcps) is one of the seventeen groups of CTLDcps discovered in mammals and composed by four members: CD93, Clec14A, CD248 and Thrombomodulin, which have shown to be important players in cancer and vascular biology. Although these proteins belong to the same family, their phylogenetic relationship has never been dissected. To resolve their evolution and characterize their protein domain composition we investigated CTLDcp genes in gnathostomes and cyclostomes and, by means of phylogenetic approaches as well as synteny analyses, we inferred an evolutionary scheme that attempts to unravel their evolution in modern vertebrates. RESULTS: Here, we evidenced the paralogy of the group XIV of CTLDcps in gnathostomes and discovered that a gene loss of CD248 and Clec14A occurred in different vertebrate groups, with CD248 being lost due to chromosome disruption in birds, while Clec14A loss in monotremes and marsupials did not involve chromosome rearrangements. Moreover, employing genome annotations of different lampreys as well as one hagfish species, we investigated the origin and evolution of modern group XIV of CTLDcps. Furthermore, we carefully retrieved and annotated gnathostome CTLDcp domains, pointed out important differences in domain composition between gnathostome classes, and assessed codon substitution rate of each domain by analyzing nonsynonymous (Ka) over synonymous (Ks) substitutions using one representative species per gnathostome order. CONCLUSIONS: CTLDcps appeared with the advent of early vertebrates after a whole genome duplication followed by a sporadic tandem duplication. These duplication events gave rise to three CTLDcps in the ancestral vertebrate that underwent further duplications caused by the independent polyploidizations that characterized the evolution of cyclostomes and gnathostomes. Importantly, our analyses of CTLDcps in gnathostomes revealed critical inter-class differences in both extracellular and intracellular domains, which might help the interpretation of experimental results and the understanding of differences between animal models.

Thyroid and endostyle development in cyclostomes provides new insights into the evolutionary history of vertebrates.

BACKGROUND: The endostyle is an epithelial exocrine gland found in non-vertebrate chordates (amphioxi and tunicates) and the larvae of modern lampreys. It is generally considered to be an evolutionary precursor of the thyroid gland of vertebrates. Transformation of the endostyle into the thyroid gland during the metamorphosis of lampreys is thus deemed to be a recapitulation of a past event in vertebrate evolution. In 1906, Stockard reported that the thyroid gland in hagfish, the sister cyclostome group of lampreys, develops through an endostyle-like primordium, strongly supporting the plesiomorphy of the lamprey endostyle. However, the findings in hagfish thyroid development were solely based on this single study, and these have not been confirmed by modern molecular, genetic, and morphological data pertaining to hagfish thyroid development over the last century. RESULTS: Here, we showed that the thyroid gland of hagfish undergoes direct development from the ventrorostral pharyngeal endoderm, where the previously described endostyle-like primordium was not found. The developmental pattern of the hagfish thyroid, including histological features and regulatory gene expression profiles, closely resembles that found in modern jawed vertebrates (gnathostomes). Meanwhile, as opposed to gnathostomes but similar to non-vertebrate chordates, lamprey and hagfish share a broad expression domain of Nkx2-1/2-4, a key regulatory gene, in the pharyngeal epithelium during early developmental stages. CONCLUSIONS: Based on the direct development of the thyroid gland both in hagfish and gnathostomes, and the shared expression profile of thyroid-related transcription factors in the cyclostomes, we challenge the plesiomorphic status of the lamprey endostyle and propose an alternative hypothesis where the lamprey endostyle could be obtained secondarily in crown lampreys.

Absence of some cytochrome P450 (CYP) and hydroxysteroid dehydrogenase (HSD) enzymes in hagfishes.

Corticosteroids are synthesized from cholesterol by steroidogenic enzyme catalysts belonging to two main families: the cytochrome p450s (CYPs) and hydroxysteroid dehydrogenases (HSDs). The action of these steroidogenic enzymes allows the genesis of the terminal active corticosteroids 11-deoxycortisol (S), 1ɑ-hydroxycorticosterone (1α-OH-B), or cortisol in different fish species. However, for Cyclostomes like hagfishes, the terminal corticosteroid is still undefined. In this study, we examined the presence or absence of CYPs and HSDs as traits in fishes to gain insight about the primary corticosteroid synthesis pathways of the hagfishes. We used published cytochrome c oxidase I (COXI) amino acid sequences to construct a phylogeny of fishes and then mapped the CYPs and HSDs as morphological traits onto the tree to predict the ancestral character states through ancestral character reconstruction (ACR). There is a clear phylogenetic signal for CYP (i.e., CYP11a1, 17, 21, and 11b) and HSD (i.e., 11-ßHSD and 3ß-HSD) derivatives of interest throughout the more derived fishes. Using trait-based ACR, we also found that hagfishes possess genes for 3ß-HSD, CYP11a1, CYP17, and CYP21. Importantly, the presence of CYP21 implies that hagfish can synthesize 11-deoxycorticosterone (11-DOC) and S. Previous research demonstrated that despite hagfish having CYP21, neither 11-DOC nor S could be detected in hagfish. This discrepancy between the presence of steroidogenic enzymes and products brings into question the expression and/or function of CYP21 in hagfishes.

Molecular identification and postprandial regulation of glucose carrier proteins in the hindgut of Pacific hagfish, Eptatretus stoutii.

Hagfish are an excellent model species in which to draw inferences on the evolution of transport systems in early vertebrates owing to their basal position in vertebrate phylogeny. Glucose is a ubiquitous cellular energy source that is transported into cells via two classes of carrier proteins: sodium-glucose-linked transporters (Sglt; Slc5a) and glucose transporters (Glut; Slc2a). Although previous pharmacological evidence has suggested the presence of both sodium-dependent and -independent transport mechanisms in the hagfish, the molecular identities were heretofore unconfirmed. We have identified and phylogenetically characterized both a Slc5a1-like and Slc2a-like gene in the Pacific hagfish (Eptatretus stoutii), the latter sharing common ancestry with other glucose-transporting isoforms of the Slc2a family. To assess the potential postprandial regulation of these glucose transporters, we examined the abundance and localization of these transporters with qPCR and immunohistochemistry alongside functional studies using radiolabeled d-[14C]glucose. The effects of glucose or insulin injection on glucose transport rate and transporter expression were also examined to determine their potential role(s) in the regulation of intestinal glucose carrier proteins. Feeding prompted an increase in glucose uptake across the hindgut at both 0.5 mM (∼84%) and 1 mM (∼183%) concentrations. Concomitant increases were observed in hindgut Slc5a1 protein expression. These effects were not observed following either of glucose or insulin injection, indicating these postprandial factors are not the driving force for transporter regulation over this timeframe. We conclude that Pacific hagfish utilize evolutionarily conserved mechanisms of glucose uptake and so represent a useful model to understand early-vertebrate evolution of glucose uptake and regulation.

Genetic and functional analysis of a Pacific hagfish opioid system.

The actions of endogenous opioids and nociceptin/orphanin FQ are mediated by four homologous G protein-coupled receptors that constitute the opioid receptor family. However, little is known about opioid systems in cyclostomes (living jawless fish) and how opioid systems might have evolved from invertebrates. Here, we leveraged de novo transcriptome and low-coverage whole-genome assembly in the Pacific hagfish (Eptatretus stoutii) to identify and characterize the first full-length coding sequence for a functional opioid receptor in a cyclostome. Additionally, we define two novel endogenous opioid precursors in this species that predict several novel opioid peptides. Bioinformatic analysis shows no closely related opioid receptor genes in invertebrates with regard either to the genomic organization or to conserved opioid receptor-specific sequences that are common in all vertebrates. Furthermore, no proteins analogous to vertebrate opioid precursors could be identified by genomic searches despite previous claims of protein or RNA-derived sequences in several invertebrate species. The presence of an expressed orthologous receptor and opioid precursors in the Pacific hagfish confirms that a functional opioid system was likely present in the common ancestor of all extant vertebrates some 550 million years ago, earlier than all previous authenticated accounts. We discuss the premise that the cyclostome and vertebrate opioid systems evolved from invertebrate systems concerned with antimicrobial defense and speculate that the high concentrations of opioid precursors in tissues such as the testes, gut, and activated immune cells are key remnants of this evolutionary role.

Evolution of variable lymphocyte receptor B antibody loci in jawless vertebrates.

Three types of variable lymphocyte receptor (VLR) genes, VLRA, VLRB, and VLRC, encode antigen recognition receptors in the extant jawless vertebrates, lampreys and hagfish. The somatically diversified repertoires of these VLRs are generated by serial stepwise copying of leucine-rich repeat (LRR) sequences into an incomplete germline VLR gene. Lymphocytes that express VLRA or VLRC are T cell-like, while VLRB-expressing cells are B cell-like. Here, we analyze the composition of the VLRB locus in different jawless vertebrates to elucidate its configuration and evolutionary modification. The incomplete germline VLRB genes of two hagfish species contain short noncoding intervening sequences, whereas germline VLRB genes in six representative lamprey species have much longer intervening sequences that exhibit notable genomic variation. Genomic clusters of potential LRR cassette donors, fragments of which are copied to complete VLRB gene assembly, are identified in Japanese lamprey and sea lamprey. In the sea lamprey, 428 LRR cassettes are located in five clusters spread over a total of 1.7 Mbp of chromosomal DNA. Preferential usage of the different donor cassettes for VLRB assemblage is characterized in our analysis, which reveals evolutionary modifications of the lamprey VLRB genes, elucidates the organization of the complex VLRB locus, and provides a comprehensive catalog of donor VLRB cassettes in sea lamprey and Japanese lamprey.

Ancient fishes and the functional evolution of the corticosteroid stress response in vertebrates.

The neuroendocrine mechanism underlying stress responses in vertebrates is hypothesized to be highly conserved and evolutionarily ancient. Indeed, elements of this mechanism, from the brain to steroidogenic tissue, are present in all vertebrate groups; yet, evidence of the function and even identity of some elements of the hypothalamus-pituitary-adrenal/interrenal (HPA/I) axis is equivocal among the most basal vertebrates. The purpose of this review is to discuss the functional evolution of the HPA/I axis in vertebrates with a focus on our understanding of this neuroendocrine mechanism in the most ancient vertebrates: the agnathan (i.e., hagfish and lamprey) and chondrichthyan fishes (i.e., sharks, rays, and chimeras). A review of the current literature presents evidence of a conserved HPA/I axis in jawed vertebrates (i.e., gnathostomes); yet, available data in jawless (i.e., agnathan) and chondrichthyan fishes are limited. Neuroendocrine regulation of corticosteroidogenesis in agnathans and chondrichthyans appears to function through similar pathways as in bony fishes and tetrapods; however, key elements have yet to be identified and the involvement of melanotropins and gonadotropin-releasing hormone in the stress axis in these ancient fishes warrants further investigation. Further, the identities of physiological glucocorticoids are uncertain in hagfishes, chondrichthyans, and even coelacanths. Resolving these and other knowledge gaps in the stress response of ancient fishes will be significant for advancing knowledge of the evolutionary origins of the vertebrate stress response.

Connexins during 500 Million Years-From Cyclostomes to Mammals.

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Novel developmental bases for the evolution of hypobranchial muscles in vertebrates.

BACKGROUND: Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic. RESULTS: Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle. CONCLUSIONS: We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Hagfish from the Cretaceous Tethys Sea and a reconciliation of the morphological-molecular conflict in early vertebrate phylogeny.

Hagfish depart so much from other fishes anatomically that they were sometimes considered not fully vertebrate. They may represent: (i) an anatomically primitive outgroup of vertebrates (the morphology-based craniate hypothesis); or (ii) an anatomically degenerate vertebrate lineage sister to lampreys (the molecular-based cyclostome hypothesis). This systematic conundrum has become a prominent case of conflict between morphology- and molecular-based phylogenies. To date, the fossil record has offered few insights to this long-branch problem or the evolutionary history of hagfish in general, because unequivocal fossil members of the group are unknown. Here, we report an unequivocal fossil hagfish from the early Late Cretaceous of Lebanon. The soft tissue anatomy includes key attributes of living hagfish: cartilages of barbels, postcranial position of branchial apparatus, and chemical traces of slime glands. This indicates that the suite of characters unique to living hagfish appeared well before Cretaceous times. This new hagfish prompted a reevaluation of morphological characters for interrelationships among jawless vertebrates. By addressing nonindependence of characters, our phylogenetic analyses recovered hagfish and lampreys in a clade of cyclostomes (congruent with the cyclostome hypothesis) using only morphological data. This new phylogeny places the fossil taxon within the hagfish crown group, and resolved other putative fossil cyclostomes to the stem of either hagfish or lamprey crown groups. These results potentially resolve the morphological-molecular conflict at the base of the Vertebrata. Thus, assessment of character nonindependence may help reconcile morphological and molecular inferences for other major discords in animal phylogeny.

Inner ear development in cyclostomes and evolution of the vertebrate semicircular canals.

Jawed vertebrates have inner ears with three semicircular canals, the presence of which has been used as a key to understanding evolutionary relationships. Ostracoderms, the jawless stem gnathostomes, had only two canals and lacked the lateral canal1-3. Lampreys, which are modern cyclostomes, are generally thought to possess two semicircular canals whereas the hagfishes-which are also cyclostomes-have only a single canal, which used to be regarded as a more primitive trait1,4. However, recent molecular and developmental analyses have strongly supported the monophyly of cyclostomes5-7, which has left the evolutionary trajectory of the vertebrate inner ear unclear8. Here we show the differentiation of the otic vesicle of the lamprey Lethenteron camtschaticum and inshore hagfish Eptatretus burgeri. This is the first time, to our knowledge, that the development of the hagfish inner ear is reported. We found that canal development in the lamprey starts with two depressions-which is reminiscent of the early developmental pattern of the inner ear in modern gnathostomes. These cyclostome otic vesicles show a pattern of expression of regulatory genes, including OTX genes, that is comparable to that of gnathosomes. Although two depressions appear in the lamprey vesicle, they subsequently fuse to form a single canal that is similar to that of hagfishes. Complete separation of the depressions results in anterior and posterior canals in gnathostomes. The single depression of the vesicle in hagfishes thus appears to be a secondarily derived trait. Furthermore, the lateral canal in crown gnathostomes was acquired secondarily-not by de novo acquisition of an OTX expression domain, but by the evolution of a developmental program downstream of the OTX genes.

Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates.

Hox genes exert fundamental roles for proper regional specification along the main rostro-caudal axis of animal embryos. They are generally expressed in restricted spatial domains according to their position in the cluster (spatial colinearity)-a feature that is conserved across bilaterians. In jawed vertebrates (gnathostomes), the position in the cluster also determines the onset of expression of Hox genes (a feature known as whole-cluster temporal colinearity (WTC)), while in invertebrates this phenomenon is displayed as a subcluster-level temporal colinearity. However, little is known about the expression profile of Hox genes in jawless vertebrates (cyclostomes); therefore, the evolutionary origin of WTC, as seen in gnathostomes, remains a mystery. Here, we show that Hox genes in cyclostomes are expressed according to WTC during development. We investigated the Hox repertoire and Hox gene expression profiles in three different species-a hagfish, a lamprey and a shark-encompassing the two major groups of vertebrates, and found that these are expressed following a whole-cluster, temporally staggered pattern, indicating that WTC has been conserved during the past 500 million years despite drastically different genome evolution and morphological outputs between jawless and jawed vertebrates.